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NOMENCLATURE 


A = first end-body 
B = second end-body 
C — center of mass of the tether system 
E — longitudinal stiffness of the tether 
F = non-gravitational forces 
g = gravitational acceleration 
Jo — moment of inertia of the tether system 
L = total tether length 

Lk — length of the &-th segment of the tether 
m a = mass of the first end- body 
= mass of the second end-body 
rrik = embedded mass k 
M — total mass of the tether system 
R = geocentric radius-vector 
Re — mean radius of the Earth 
Rc = geocentric radius of the center of mass 
$ — arclength along the unstretched tether 
T — tether tension 
T = tether tension vector 
i = time 

7 = tether elongation 
p = tether mass per unit length 
T = unit vector along the tether line 
T( = direction of the imaginary straight tether line 
1? = rotational angular rate 
( ’ ) — differentiation with respect to time 
( 1 ) — differentiation with respect to the arclength 
(a,b) — scalar product of vectors a and b 
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1 . EQUATIONS OF MOTION IN NEWTONIAN FORM 


In the first part of this study [1], we considered the dynamics of a momentum 
exchange system consisting of two end-bodies and a number of power stations con- 
nected with tether segments. The end-bodies A and B and the power stations k 
were modeled as point masses ni A , m B , and mi, respectively. The power stations 
are connected with tether segments of lengths L^, as shown in Fig. 1. The tether 
segment Z* connects masses m*. and rrik+i. 


L« L, ... L v 


m A mi m 2 ... m k m k+ i m B 


Fig* 1. Structure of the momentum exchange tether system. 

All point masses and the masses of the tether segments arc assumed to be 
constant. Tether mass per unit length can vary along the tether. 

Positions of the tether elements with respect to a non-rotating geocentric 
reference frame OX Y Z are defined by the geocentric radius R as a function of the 
arclength .s measured along the unstretched tether from .4 to /?, and time t , 

R = R (s, t) y 

Positions of the end masses and embedded masses are 

R a = R(s > i,t) > Rb = R(sjj,<), Rfc = R 
The tension vector T of a perfectly flexible tether is tangent to the tether line, 

T = Tr, , " 7 = |R'I, ( 1 ) 

7 

where T is a unit vector tangent to the tether line, prime denotes differentiation 
with respect to the arclength s, and 7 is the local tether elongation. 

The tether tension T can be expressed as a function of the elongation 7, 
elongation rate 7, temperature 0, and other factors, 

T — T(s, t, 7, 7, < 9 ,...). ( 2 ) 
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k 


A 


Fig. 2. Positions of the tether system elements. 

Equations of motion of the tether system derived in [1] include ordinary and 
partial differential equations. The motion of the end masses and embedded masses 
is described by the ordinary differential equations 

m A = T 4 + m A g 4 4- F a 

m B R s = -T w 4- m B g B 4- F B (3) 

m k R-k = T fc+ - Tjt_ 4- m k g* 4- F* 

where dots denote differentiation with respect to time, g A , g B , and g fc are the 
gravity accelerations at points A, B, and k, respectively, while F 4 , F B , and F*. 
are non-gravitational forces acting on the end masses and embedded masses. The 
tether tension vectors are taken at the following points: at point A, T B at point 

B, T*_ at point k of segment I*_i, and T* + at point k of segment L k . 

The motion of the tether is described by the partial differential equation 

pR = T +pg + F (4) 

where dots denote differentiation with respect to time t , and primes denote differ- 
entiation with respect to the arclength s, p is the tether mass per unit length, and 
T is the tether tension. 

In the general case, the tether mass per unit length p can vary along the 
tether, but we will concentrate on a more practical case when p is constant along 
each tether segment connecting the end masses, 

P = Pk, k = l,...,K. (5) 
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2. EQUATIONS OF MOTION IN MINAKOV’S FORM 


Equations of motion in Minakov’s form, introduced into space tether dynamics 
in [2], are obtained as follows. First, we differentiate equation (4) with respect to 
the arclength s, 

pR' = T" +pg' + F'. 

Keeping in mind that R' = 7 x and T = Tx, we find that 


P (-yT + 27T + 7x) = T"x + 2 TV + Tt" + p g' + F'. 


( 6 ) 


By definition, the unit tangent vector x satisfies the condition 


The first differentiation with respect to the arclength and time yields two orthogo- 
nality conditions, 

(t,t) = 0, (t,x') = 0, (8) 

while the second differentiation yields two kinematic relations, 

(x, x) = -(x, x), (x, x") = -(x', x'). (9) 


Multiplying equation (6) by vector x (scalar product) and taking into account 
relations (7)-(9), we obtain a scalar equation of the second order with respect to 
tension, 

P7 - P7 (t, t) - T" - T(x',x') + p(g',x) + (F',x). (10) 

From this equation, we express the second derivative of tension as 


T " = +P1 ~ Pli*,*) -p(g',x) - (F',x). (11) 

Substituting expression (11) into equation (6), we arrive at the following equation 
of the second order with respect to the unit tangent vector X, 

p [ 7 T + 27 X + 7 t (t, x)] = T [x" + (x' , x 1 ) x] + 2 TV + pg' + F' - (pg' + F' , x) x. ( 12 ) 


This equation describes transverse motion of the tether and it has an obvious 
wave structure with a transverse wave velocity equal to 


Equation (10) describes longitudinal motion of the tether and it also has a 
wave structure. To reveal this structure, we consider a case of linear elasticity, 

T = T 0 + E( 7-7,), (14) 

where T 0 is a reference tension, E is the tether longitudinal stiffness, and 7, is 
he equilibrium elongation under the reference tension T 0 . According to (14) and 
keeping in mind that 7* may vary with time and temperature, we have 

- f „ 

7 “ E + 7 *' 

After substituting this relation into equation (10), we can clearly see that this 
equation has a wave structure with a longitudinal wave velocity equal to 


(15) 



The longitudinal wave velocity (15) is much higher than the transverse wave 
velocity (13), and therefore, the system of equations (10) and (12) is stifT in the sense 
that equation (12) describes high frequency longitudinal oscillations, while equation 
(10) describes transverse oscillations of much lower frequencies. The advantage of 
Minakov’s form of the equations lies in the separation of the high and low frequency 
wave structures in the tether dynamics equations. 

In the Newtonian field, 


the gravity gradient along the tether is calculated as 


g = - 




R 3 


r_Bi 

R 3 


+ ifHR>R') = ^[3e*(e*,T)- T ], 


(16) 


where fi B is the gravitational constant of the Earth, e K = R/R is a unit vector along 
the geocentric radius vector R, and R' = 7 T. The corresponding gravitational term 
in the equation of longitudinal motion (10) is reduced to 


(S\ r )= ~^f-l3(e R ,T) 2 -1], ( 17 ) 

while the respective gravitational term in the equation of transverse motion (12) 
takes the form 


S' ~ (s',t)t= ^ 3(e H ,T)[e R - (e*,T)T]. (18) 
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3. BOUNDARY CONDITIONS 


The equations of tether motion must be complemented with boundary condi- 
tions at the ends of each tether segment. To derive these boundary conditions, we 
express the accelerations of the end masses A and B and the embedded masses* m k 
from equation (3) 


n _ T* + F 4 

— _ b g4, 


Rb = 


m A 

— T g + F B 

m B 


+ gs, 


T) T *+ - T*_ -|- 
R ‘ = +g ‘> 


and from equation (4) at the ends of the tether segments, 


‘-Ft 1 ) 


+ g4, 


+ g B, 


+ g * 




+ gfc* 


This yields the following boundary conditions 


T f |a = Q,t, T'| b = Q b , 

T'| fe _ = Qt_, T'|*+ = Q k+> 


(19) 


where 


Q* = ~~ (T4 + F,*) - F| a , 

m A 

Qb = —(-T b +F b )-F\ b , 

TTlfi 

Q*- = ^ (T* + - T*_ + Ft) - F|*_, 

m t 

Q*+ = — (T* + - Tt_ + Ft) - F|*+. 

mt 


( 20 ) 


Note that the gravitational accelerations canceled out and are not included 
in the boundary conditions. 



Keeping in rnind that T' = T't + Tt\ according to (1), we multiply each of 
the relations {19} by the corresponding unit tangent vector t (scalar product) and 
arrive at the following conditions 

r'u = (Q A ,T < ), ru = (Q B , Ti ), 

r\ u . = (q».,t*-), tv - (Q 4 + ,T i+ ), (21) 

Relations (21) serve as boundary conditions for equation (10) describing longitudi- 
nal motion of the tether. 

Substituting relations (21) into (19), we find boundary conditions for equation 
(12) describing transverse motion of the tether, 


T '\a = ~ [Qa - (Qa,t a )t a ], 

j A 

T 'U = jT [Q« - 

T 'i*- = [Qfc- - (QA.,T fc _)r fc _], 


Ah+ = 


jrHQ*+ - (Qjb4,Tfc + )T fc+ ]. 


( 22 ) 


The boundary conditions are simplified when the non -gravitational forces F 
are neglected. The boundary conditions for the t ransverse motion of the tether are 
reduced to 

Aa = o, 

As = 0 , 

Ah- = ^ (T fe+ - **T fe _), (23) 

T 'l‘+ - P ir -*»-), 

Wfe -lk+ 

where 

= (Tfe+,Tjfc_), 


while the boundary conditions for the longitudinal motion of the tether take the 
form 


T\ a = — 1 

m A 

^ 




T'U = -^-T. 


m B 




m (*kTk+ - T fe _), T% + = 

mj. 


(24) 
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4. MOTION OF THE CENTER OF MASS 


The equations of tether motion derived in the previous sections describe evo- 
lution of the tether orientation and shape. They must be solved together with one 
of the ordinary differential equations (3) describing the orbital motion of one of the 
end masses or embedded masses, or with the equation of motion of the center of 
mass 

= (®o + 4>o), (25) 

where R c is the geocentric radius vector of the center of mass, M is the total mass 
of the tether system, G 0 is the sum of the gravitational forces 

rB 

G 0 = m A g„ + m B g B + ^ m/bgfc + / pg ds , (26) 

k 

and <|) 0 is the sum of the non-gravitational forces acting on the tether system, 

r B 

<I> 0 = F 4 F b + Y, F* + / F ds. (27) 

k 

As shown in the first part of this report [1], the sum of the gravitational forces 
can be represented as a series 


oo 

Go = Mg“ + (28) 

n=2 

where g° is the gravitational acceleration at the center of mass, g" are derivatives 
of the gravitational acceleration along the “mean” tether line R = R c + (.? - Sc ) Tl 
drawn through the center of mass, 


n= 1 

” c n! ds n 


and quantities I n represent the high order moments 

rB 

I n — fn A (sx — «c) n + (s B — Sc) n + ^ m* ( Sk — $c) n + / p(s — sc) n ds. 

t Ja 
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The direction of the “mean’ tether line T t is defined as follows. As shown in 
the first part of this report [1], the tether shape can be represented by the series 


R(.s,t) _ R c (t) + (s — «c)T|(f) + ^ q>t(0^n(s), (29) 

»= 2 

where q n (t) are generalized coordinates and U n (s) are the eigenforms of tether 
oscillations. Differentiating this equation with respect to the arclength, and keeping 
in mind that R/ = ^x, we derive 


T » =TT- J]q„(<)^n(*)- (30) 


The orthogonality conditions for the eigenforms can be expressed as 

J A PVlV' j ds = 0, t ^ j, (31) 

where P = P($) represents a normalized tension profile, defined as the solution to 
the following boundary problem [1], 


P ' = -/>(«- s c ), 

Pa = ~m A (s A - s c ), 

Pb = m B (s B - s c ), 

Pk+ = Pk- - m k (sk - s c ). 


(32) 


The first eigenform describes the rigid rotation £/, = s-s c , for which U' = 1. 
Multiplying equation (30) by PU[ = P , integrating over the entire tether length, 
and applying the orthogonality conditions (31), we find that 

B 

Pyrds. (33) 

i 

The integral of P , according to [1], is equal to the moment of inertia of the unde- 
formed tether system about the center of mass 



/ Pds = \\U l \\ = J c = 

J A 
f B 

m A (sa ~ s c ) 2 + m B (s B - Sc) 2 + ^ m* (s* — s c ) 2 + I p(s — s c ) 2 ds, 

k Ja 


(34) 


46 


and therefore, relation (33) can be rewritten as 


Ti 


Jc 



Pyrds. 


(35) 


Formula (35) defines the “mean” tether direction for any given tether shape under 
the assumptions of the current formulation. 

Knowing T, , we can calculate derivatives (26). As shown in [1], the Newtonian 
term of the gravitational field gives 


8c = 


PE 7," 

Rl +n 


[ e fl P n+l( x ) + e > ^n( a: )] 


(36) 


where c R = R C /|R C | is the unit vector of the geocentric direction to the center of 
mass, 7 ! = |Ti|, e x = Ti/h, P n {x) are Legendre polynomials, and x = ~(e fi ,ei). 

To close the system of equations, we need to find the geocentric positions of 
all elements of the tether system. For a given tether shape, we know the relative 
positions r = R — R* of all elements, and from the definition of the center of mass, 
we can find its relative position 

r c — + 'y + I prdsj. (37) 


Using this information, we can calculate the geocentric position of the end 
point R a — R c — t c and the geocentric positions of all other points of the tether 
system R = R^ + r. 


5. QUASI-STATIC TENSION 


The longitudinal oscillations of the tether system are directly affected by the 
internal friction in the tether, while the transverse oscillations are are affected 
indirectly, through a weak coupling with the longitudinal oscillations. 

As discussed in [2], one way to model internal friction in the tether is to add 
a term px 7, where \ is an effective damping coefficient, to the left side of the equa- 
tion of longitudinal motion (10). It is shown in [2] that the energy dissipation in 
the tether can rapidly damp out high frequency longitudinal oscillations, and after 
a short period of t im e the longitudinal motion becomes quasi-static. Tension vari- 
ations in the quasi-static motion are induced by relatively slow transverse motions 
and can be described by equation (11) without the term p 7 , 

r - T(V, t') - - p (g'> t) - (F». ( 38 ) 
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At any given moment f, this equation can be treated as an ordinary differential 
equation with respect to the tether tension T. Along with the boundary conditions 
i 21 ?, it defines the tether tension profile as a function of the dynamic state of the 

tether, including the tether shape, elongation and the angular rate distribution 
along the tether. 

Similar to (24), the boundary conditions ( 21 ) can be presented as 


T'U = ^T i + / A , 

m A 

T'\ B = -£*- Ta +f B , 

Til 

T> ^+ = ( Tk + - + fk+, 

Tflfc 


(39) 


where 


ffl A 

* = r.,T»)-(r | „T,). 

/*-= ~ (F»,t,-)-(FU-,t*_), 

At = (Fi>1»+) -(F|t+,T, + ). 


As we see in relation ( 16 ) , the gravity gradient g' calculated along the tether 
is proportional to the elongation 7 , g' - 7 g'. This is true for any gravitational 
field, not only Newtonian. Also, it has been shown in [ 2 ] that the aerodynamic 
and Ampere forces, as well as solar radiation pressure acting on a tether element 
are proportional to the elongation, F = 7 F 0 . Formally, the gradient of these forces 
F' will include quadratic terms in 7 because of the variation of the environmental 
parameters along the tether, however, these terms can be linearized with respect 
to small variations of elongation. 

Under these assumptions, if we use the linear elasticity relation (14) between 
the elongation and tension, then we will have 


7 = 7 , + (T - T 0 )/E, (40) 

and equation (38) will be linear with respect to T . The boundary conditions ( 39 ) 
will also be linear with respect to T. 

The general solution of the linear boundary problem (38)— f 39) at any given 
moment t can be represented as a sum of two solutions, 

^ (-SO — ^o( s jt) + CiTi(s,t). (41) 
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Solution T 0 is obtained as follows. We set T 0A to zero and determine T' from 
the first equation of the system (39), 

T 0 A = 0, Tq A = f A . (42) 

1 hen, we integrate equation (38) with the initial conditions (42) over the first tether 
segment, and find and T' k _ at the end of the segment. Now, we nse the third 

equation of the system (39) to find the tension at the beginning of the next tether 
segment 

Ti+ = 7T k [ r ‘- + jt (r *- “*-)]• («) 

We then substitute quantity (43) into the last equation of the system (39) to find 
the derivative at the beginning of the next segment. 

This process is repeated until we arrive at the end B and determine the values 
of T cb and T‘ IB . Generally, these values will not satisfy the second equation of the 
system (39). We need solution T[ to match the boundary condition at the end B. 

Solution 7\ is calculated as follows. We introduce a new system of equations 

(38’ j (39') which retains only terms linear in T from the original equations (38)- 

(39). We set T lA to some reference tension , 

T,a=T„ ( 44 ) 

find T; a from the first equation of the homogeneous system (39'), and integrate the 
homogeneous equation (38') over the first tether segment. We then use the third 
and the last equation of system (39') to find 2*+ and T‘ k+ at the beginning of the 
next segment, and continue integration until we reach the end B. 

Now, we can substitute the general form of solution (41) into the second 
equation ol (39) and find the coefficient Cj from the following equation 

+ m. *B Tr3B + Cl ( T ' b + ^ TiB ) = /b ’ ( 4S ) 

One of the major benefits of this approach is that we reduced the initially stiff 
system of equations (10), (12) to a non-stiff system (12), (38) by abstracting from 
the fast transient processes in longitudinal oscillations. This does not mean that 
we neglected tether elongation, we only filtered out its high frequency components. 


6. SIMULATION APPROACH 


For a numerical solution, equations (12), (22), (38), (39) can be discretized 
with respect to the arclength s. This will reduce the partial differential equation 
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(12) to a set of ordinary differential equations at the discretization nodes. These 
equations are complemented with the ordinary differential equation (25) describing 
the motion of the center of mass. 

The number of discretization nodes is determined by the required accuracy. 
To get an idea of a typical number of nodes, let us consider a very simple model of 
a siring clamped at the ends A and H. Its transverse oscillations are decribed by 
the string equation 

pil — Tu'\ 

where p is the string mass per unit length, T is the tension, assumed constant 
along the string, u = u(s,f) is the transverse displacement, dots denote differenti- 
ation with respect to time i, and primes denote differentiation with respect to the 
ardength s. 

The boundary conditions, are, obviously, 

Ua — Ub = 0 , 


A discretized string equation looks as follows 

2 u k -\- 1 ~ + Ufa— i 


u = v: 


as* 


where v t - y/Tjp is the transverse wave velocity, ti* is the transverse displacement 
at the node fc, as is the distance between the neighboring nodes, 


AS = 


N * 


N is the number of nodes, and L is the length of the string. 

Partial solutions representing the natural modes of oscillations are found in 
the form 

u k - sin(J?„f) sin(Tm^), n = 1,2, . . . 

where O n are the eigenfrcquencies, and — k /N, k = 0, 1, , N . 

Substituting this form of solution into the discretized string equation, we find 
that the eigenfrcquencies of the discretized problem are equal to 

n ojir Vt • f‘ Kn \ v * (-, 

L \2N J L y 24 N 2 J 

The exact eigenefrcquency obtained from the original partial differential equa- 
tion is equal to 

rt* Vt 

il n — irn — * 

L 


50 



The relative difference is 


^n. ~ ^ 1 /irn\ z 

fl* n ~ 24 \JV"/ * 

Because of the difference in eigene frequencies, the phase difference between 
the simulated and exact solutions will accumulate with time as ay? = (Q n - £2' n )t. 
If we need to simulate the excitation of the n-th mode with a relative accuracy e T1 

over a period of time af, then the difference between the eigenefrequencies must be 
limited as 

IA. - All < % 

At 

To satisfy this relation, the number of nodes must be at least 


N > 

nj 



For a typical momentum exchange system with L = 90 km and v t = 1 krn/s, 
it takes at least 125 nodes to maintain an accuracy of 1% for the first mode of 
transverse oscillations over one 3-hour orbit. 

If an explicit integration scheme is used, then the time step will be limited by 
the condition of stability of the numerical solution. In general, this restriction has 
a form of 


M < kt min ( — — 1, 
* 1 ®ti J ' 


where k t is a numeric coefficient depending on the implementation, as; are the 
distances between the nodes, and vn are the maximum transverse wave velocities 
(13) between the corresponding nodes. 

By comparison, in a lumped mass model, the time step is limited by 


A f> ^ 1?^ 



where k r is a numeric coefficient depending on the implementation, and v ei are the 
maximum longitudinal wave velocities (13) between the corresponding nodes. 

The longitudinal wave velocity is much higher than the transverse wave ve- 
locity, and therefore the time step will be much smaller. This is why simulations 
based on the Minakov’s form of the equations of motion are much faster than those 
based on lumped mass models. 

Two numerical models were implemented, one using Minakov’s form of the 
equations, and the other using a lumped mass model. With a sufficient number of 
nodes, as explained above, the results of the simulations were consistent between 
the models within the required accuracy. As expected, the Minakov’s type model 
was by an order of magnitude faster than the lumped mass model. 
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7. RESONANT EXCITATION 


It has been noted in simulations that the transverse oscillations of a typical 
momentum exchange system go through cycles of a long-term amplitude modula- 
tion, when the amplitude first builds up during a number of consecutive perigee 
passages, and then decreases again. 

Fig. 3 shows a history of the transverse amplitude variations over the course 
of twelve orbits. The amplitude profile is shown in light gray, while the perigee 
passages are marked with black clusters underneath. 



Fig. 3. Resonant excitation of the transverse oscillations. 

We see that even without the initial excitation the amplitude of transverse 
oscillations steadily grows during five perigee passages until it reaches a relatively 
high level of over 120 m T stays at this level for another perigee passage* and then 
steadily decreases during the next five perigee passages. This cycle repeats with 
minor variations. 

An in-depth analysis of this phenomenon based on the modal formulation de- 
veloped in the first part of this report reveals that the mechanism of this excitation 
is resonant in nature and is caused by a sharp tuning between the natural frequency 
of the main mode of transverse oscillations U 2 and the frequency of the gravity gra- 
dient variation, which cycles twice per every rotation of the tether system about 
its center of mass. 

This timing is persistent in the sense that the parameters of the tether system 
must be changed radically to get away from the resonance. In particular* changing 
the spin rate has very little effect because the transverse frequencies are proportional 
to the spin rate, as shown in [1]. 
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The modal formulation of Part f accurately describes the motion of the tether 
system in general and the resonant excitation of transverse oscillations in particular 
when the term Q* in the main modal equation (37) derived in Part I is used in its 
general form (38). The simplifications of the term Q n demonstrated for a special 
(.ase o an ideally tapered tether should remain restricted to this special case and 

should not be applied in case when the MXER tether is not carrying the full load 
lor which it was designed* 

tor practical purposes, the modal solution requires very few modes and is 
computed much faster than the numeric solutions based on the Minakov’s formu- 
latiGii, not to mention lumped mass models. 


8* SENSITIVITY TO NON-GRAVITATIONAL PERTURBATIONS 
8.1* Method of Computation. 

To evaluate the sensitivity of the motion of the momentum exchange system to non- 

gravitational perturbations, we use the modal decomposition method described in 
Part I of this study [1]. 

As shown in [1], the effect of non-gravitational forces on the excitation of the 
eigen form V n is determined by a generalized force 

cp n = F A U nA + F B U nB + U nk + f FU n ds. (46) 

k Ja 

It is computationally expensive to calculate the integral of F U n over the entire 
tether length on each step. To make this calculation more efficient, we will use the 
following quadratic approximation on each tether segment s k < s < s k+1 

= (2f -e)F + (1 - 4£ 2 )F(s fe+x/2 ,f) + (f + 2f a )F (s k+u t), (47) 

where £ (s f 2 )/(■**-(- 1 — an d Sfe-1-1/2 = (•*& T a *+i)/2 is the middle of 

the segment, so that -1/2 < £ < 1/2, 

Using approximation (47), we have 


f FUnds = y 

JA t 


yZ Qnk, 


where 

rk+l 

Qnfc=/ F U n ds = 

J k 

i 2 *nk - 4iiJF(sfc,f) + {I° k - 4Il k )F(s h+1/2> t) + (/^ + 2/^ A .)F(s fc+ i,t), 
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and 


T° 

*nk 


rk-n 

Jk 


U n ds, 


rM-i 

= L 


f- J ' n f 


1 2 


nk 


fk+l 

Jk 


U n ( 2 ds. 


SSSAZZZZZ' - ~ d — - =s 

With tjpical parameters, the accuracy of this presentation of the generalized 

forces meets ° r as «**«-«*. ** *££££ 
dynamics of the momentum exchange system. 

8.2. Aerodynamic Forces. 

matedT" b 121 ’ tl,e “ rodynMr,ic for " “‘“g on » ‘other segment can be appro*!- 


P = -PtldtJ 


|T xv.|( (l + |) v a - 1 £Vt j +*/(!- e)Jv w j(v a - v r ) , 


(49) 


where d t is the tether diameter, p a is the air density, v c is the tether velocity relative 

a °, hc "5 Vr " ( v - i y) T 18 the component of the velocity v ft along the tether line 
and e and i/ are small parameters, 0 < e, u < 0 1 

The aerodynamic forces acting on the end masses and the embedded masses 
Ere calculated as for regular satellites [4]* 

, |’ or ^ lloytether [5] consisting of separate strands kept apart from each other 
formula (49) must be applied to all individual strands. One must also take into 
account the overshadowing of the strands. In general, a Hoytethcr will have a larger 
total exposed area of the tether and a higher air drag. 

With typical system parameters and Hoytethers having between 16 and 24 
strands, it has been determined in numerical simulations that a 5% variation in the 
aerodynamics forces due to air density variations or uncertainties of the aerody- 
namic parameters of the tether will cause a 2 m shift of the tether tip A after one 
orbit . 


8.3. Ampere Forces. 

Tlie Ampere force acting on a tether segment can be represented 

F = Ijt x B 


as 


(50) 


vector 1 ' 5 th<! elCC ‘ riC " lrient “ ‘ hC tether ’ "" d B “ the ^““gnetic induction 
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Normally, Ampere forces should not be applied during a preparation for a 
rendezvous. However, to evaluate the sensitivity to the variations of the Ampere 
forces, we will consider a purely hypothetical situation when there is an electrical 
“leak” in the system, which results in a very small alternate electric current of a 
1 mA amplitude flowing in the conductive segments of the tether in the direction 
of the EMF and modulated proportionally to the EMF. It was also assumed for 
simplicity that the “leak” current shuts off at altitudes above 2000 km for the lack 
of electrons. 

It has been observed in simulations that this electric “leak” during the perigee 
passage causes a 3 m shift of the tether tip A after one orbit. 


8.4. Solar Radiation Pressure. 


The force of the solar radiation pressure acting on a round tether can be approxi- 
mated as [2] 


F = -p t d t 7 |t x e,| 


V, 

4 . , 1 

i( i+ j) 

e, - - *(e,,T)T 


(51) 


where e, is a unit vector of the direction to the Sun, p 3 is the solar radiation 
pressure, and x is the reflection factor. 

The solar radiation pressure forces acting on the end masses and the embedded 
masses are calculated as for regular satellites [4] . 

For a Hoytether [5] consisting of separate strands kept apart from each other, 
formula (51) must be applied to all individual strands, taking into account their 
overshadowing. In general, a Hoytether will have a larger total exposed area of the 
tether. 

With typical system parameters and Hoy tethers having between 16 and 24 
strands, it has been determined in numerical simulations that a 5% variation in the 
forces of solar radiation pressure due to uncertainties of the reflective and geometric 
parameters of the tether will cause a 1 to 4 m shift of the tether tip A after one 
orbit. 


8.5. Thermal Expansion. 

Thermal expansion changes the length of the tether segments and thus changes 
the moment of inertia of the momentum exchange system. The variations of the 
moment of inertia result in variations of the tether spin rate. 

The moment of inertia is proportional to the second power of the average 
elongation 7, , and in a simple case of a free rotation, when the angular momentum 
is conserved, the angular velocity Q satisfies the relation 

7?^ = 7o^°* (52) 
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When the average elongation 7, varies with the temperature 0 as 


7i=7o + ot t a 0 , (53 ) 

where a t is the thermal expansion coefficient, then the angular rate will vary as 


7o 


7o 


If the average deviation of the actual temperature from the modeled temper- 
ature is a<9 ot , during a time af, then the angle of the tether in-plane orientation 
will deviate by ai? = and the transverse displacement of the tether end ,4 

will be estimated as 

„ r 2a ( a(9 at , 

Al) L AC fti — - do At L ac , 

where L AC is the distance from the end A to the center of mass C. 

According to [3], Zylon, a candidate tether material, has a negative thermal 
expansion coefficient of a, = -6x 10“ 6 1/K. With a typical spin period of 6.3 min 
and a distance from the end to the mass center of 75 km, we find that after one 
orbital period of M =: 3 hours, the tether end will drift away from its anticipated 
position by about 163 m per each Kelvin of the average temperature deviation. 
This sensitivity is attributed to the high spin rate Q 0 . 

8 .6. Creep. 

The average elongation 71 will gradually increase with time due to creep, 

7i=7o + « e ^, (54) 

where a c is the creep coefficient. In the approximation (52), the average spin rate 
will be dropping as 

/o 2^7, 2 a c At „ 

« - — — o Q = — — n 0 . 

7o 7o 

and the deviation of the angle of the tether in-plane orientation will be growing as 

ad fe - — — 

7o 

The transverse displacement of the tether end A will build up with time as 

< 1 ^ L ac ^ L A c } 

7o 

where l A c is the distance from the end A to the center of mass C. 
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As reported in [3], the creep coefficient a c of Zylon-HM under high stress 
is on the order of 10^ fi 1/hour, For a 100- km tether, this will result in tether 
elongation at a rate of 2*4 m/day* With a typical spin period of 6.3 min and a 
distance from the end to the mass center of 75 km, we find that after one orbital 
period of aI = 3 hours, the tether end will drift away from its anticipated position 
by about 41m. For Zylon-AS, the creep coefficient and the tether end displacement 
are 3 times higher* As with thermal expansion, this sensitivity is attributed to the 
high spin rate 

8,7* Mass Loss. 

There are several processes that will contribute to the variation of the mass of the 
momentum exchange system: 

(1) outgassing 

(2) sublimation 

(3) micrometeorite damage 

(4) hollow cathode emission 

(5) molecular deposition (mass addition) 

While there is no reliable data on most of them, it has been estimated that a typical 
momentum exchange system may lose between 5 and 150 grams of its mass per day. 
As a result, of this process, the center of mass of the tether system will be slowly 
shifting along the tether. 

Let us consider the motion of a system with a straight, quasi-rigid tether, 
whose elements m* are positioned at points along the line 


Ri = Re? + (■Si — (55) 

where x is a unit vector, and subscript C refers to the center of mass. By the 
definition of the center of mass, 

5^ mi ( si - s c ) = 0. (56) 

i 

According to (56), when the masses of the elements change, the center of mass 
slides along the tether at a rate of 

*c = m (57) 

i 

where M is the total mass of the tether system, 

M = ^ rrti. 

I 
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Differentiating (55) with respect to time, we obtain 


R, — Re — &c t — 2 Sc t •+• (it — ic) t. (58) 

Multiplying (58) by m i and performing summation, we find that 

^miRi = M(R C - ict — 2i c t), (59) 

ft 

where the terms with T canceled out because of the definition of the center of mass 
(56). 

For each element of the system, 

m»Ri = T, + F, + W it (60) 

where T,, F,-, and W,- are the sums of all internal, external, and reactive forces 
acting on this element. Substituting equations (60) for the individual elements into 
(59), we arrive at the following equation of motion of the center of mass 

R-c = ic T + 2 i c T + — ^^(F< + Wj). (61) 

ft 

All internal forces T* cancel out after summation, as dictated by the Newton^s third 
law. The terms with the derivatives of s c describe the perturbation of the orbital 
motion because of the relocation of the center of mass* 

Multiplying (58) by T (vector product) and performing summation, we find 

that 

Jc ft = 5>i -«c)tx(F, + W ( ), (62) 

ft 

where J c is the moment of inertia about the center of mass 

J c = - sc) 2 , 

t 

and O is the angular velocity of the tether rotation 

Q = T X T. 

It is interesting to note that even though the moment of inertia J c is changing 
because of the mass loss, it does not directly result in the variation of the angular 
velocity, unlike the case of creep and thermal expansion. 

While the hollow cathodes are expected to create a small torque of reactive 
forces Wj, it is not obvious that the other processes involved in the mass exchange 
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with the environment will be anisotropic enough to create any noticeable torque of 
the reactive forces. 

Assuming that the reactive torques are negligible, the gravitational field is 
Newtonian, and the non-gravitational forces are small, the in-plane perturbation of 
the orbit of the center of mass can be described, according to (61), as 


Sx — 2u) Sy — tb Sy — 


Sy -I- 2u ; Sx + cj Sx — 



2 /*g \ 

R%) 

R%) 


Sx = s c cos 0 — 2 s c f? sin t?, 
Sy = $c sin i? 4* 2sc cos i?, 


(63) 


where Sx and Sy are the vertical and horizontal displacements of the, center of mass, 
u) is the orbital angular rate, f? is the spin rate, and t? is the angle between the 
tether and the local vertical. For a fast spinning tether, J? u;, we have i? % f2t, 
and the order of magnitude of the perturbation is estimated as 


Sx 2 + Sy 2 



These displacements are expected to be very small (less than a millimeter) because 
of the fast spin and slow shift of the center of mass. 


9. ESTIMATION AND CONTROL REQUIREMENTS 


As shown in the previous section, the position of the tether tip in the momen- 
tum exchange system is very sensitive to variations of the system and environmental 
parameters, when considered from the standpoint of the rendezvous precision re- 
quirements. 

The variations producing the most impact on the tether tip positioning are 
related to the lengths of the tether segments, as we have seen in the case of thermal 
expansion and creep. 

If the distances between the neighboring tethered modules could be measured 
with an accuracy of a few millimeters, then we could use this information to sig- 
nificantly improve the estimate of the current state of the momentum exchange 
system. 

Every dynamic process in the momentum exchange system, one way or an- 
other, leaves its signature in the longitudinal motion of the tether. Using a detailed 
dynamic model, we may be able translate a continuous stream of precise distance 
measurements, combined with other data, into a reasonably accurate estimation of 
the tether system state. 


Furthermore, using the information about the system state, we could com- 
pensate for some unexpected (or unmodeled) deviations from the projected path. 
In particular, deviations of the rotation about the center of mass could be corrected 
by slightly varying the tether length(s) or changing the mass distribution of the end 
bodies or power stations. 

Given extremely short rendezvous windows with very fast relative motions, 
and various inherent uncertainties of the environment, it is very clear that elaborate 
estimation and control algorithms must be developed to make successful rendezvous 
with the momentum exchange system possible. 


10. CONCLUSIONS 


The task of prediction of the motion of a momentum exchange system with 
a 9-digit accuracy is extremely challenging, and it is relentlessly testing our ability 
to gain fundamental insights into the nature of the tether dynamics in this system. 

It has been shown that the modal decomposition approach developed in the 
first part of this study is a very powerful and precise tool for the simulation of the 
dynamics of momentum exchange tethers. 

However, to get a prediction with the required accuracy using this tool, one 
must have precise inputs, including the initial state, system parameters, and the 
environmental models. In simulations, the tether tip positioning has been observed 
to be quite sensitive to even small variations of the parameters involved in the 
calculations. 

It is therefore imperative that precise estimation and control algorithms be 
developed to compensate for the inevitable uncertainties and support high precision 
rendezvous. 
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